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Abstract 

Background: Aluminium (Al) toxicity is considered to be one of the major constraints affecting crop productivity 
on acid soils. Being a trait governed by multiple genes, the identification and characterization of novel transcription 
factors (TFs) regulating the expression of entire response networks is a very promising approach. Therefore, the aim 
of the present study was to clone, localize, and characterize the TaSTOPl gene, which belongs to the zinc finger 
family (Cys2His2 type) transcription factor, at molecular level in bread wheat. 

Results: TaSTOPl loci were cloned and localized on the long arm of homoeologous group 3 chromosomes 
[3AL (TaSTOPl -A), 3BL (TaSTOPl -B) and 3DL (TaSTOPl-D)] in bread wheat. TaSTOPl showed four potential zinc finger 
domains and the homoeologue TaSTOPl -A exhibited transactivation activity in yeast. Expression profiling of TaSTOPl 
transcripts identified the predominance of homoeologue TaSTOPl-A followed by TaSTOPl-D over TaSTOPl-B in root 
and only predominance of TaSTOPl -A in shoot tissues of two diverse bread wheat genotypes. Al and proton (H + ) stress 
appeared to slightly modulate the transcript of TaSTOPl homoeologues expression in both genotypes of bread wheat. 

Conclusions: Physical localization of TaSTOPl results indicated the presence of a single copy of TaSTOPl on 
homoeologous group 3 chromosomes in bread wheat. The three homoeologues of TaSTOPl have similar genomic 
structures, but showed biased transcript expression and different response to Al and proton (H + ) toxicity. These results 
indicate that TaSTOPl homoeologues may differentially contribute under Al or proton (H + ) toxicity in bread wheat. 
Moreover, it seems that TaSTOPl-A transactivation potential is constitutive and may not depend on the presence/ 
absence of Al at least in yeast. Finally, the localization of TaSTOPl on long arm of homoeologous group 3 
chromosomes and the previously reported major loci associated with Al resistance at chromosome 3BL, through QTL 
and genome wide association mapping studies suggests that TaSTOPl could be a potential candidate gene for 
genomic assisted breeding for Al tolerance in bread wheat. 

Keywords: Aluminium, TaSTOPl, Triticum aestivum L, in situ hybridization, Transcription factor, Transactivation, 
Homoeologue, pH 



^BMC 

Plant Biology 



* Correspondence: lishabaihe@hotmail.com; plopes@utad.pt 
Centra de Genomica e Biotecnologia, Institute de Biotecnologia e 
Bioengenharia (CGB/IBB), Universidade de Tras-os-Montes e Alto Douro 
(UTAD), P.O. Box 1013, 5001-801 Vila Real, Portugal 

2 Departamento de Genetica, Facultad de Biologia, Universidad Complutense 
de Madrid (UCM), Madrid 28040, Spain 

Full list of author information is available at the end of the article 

O© 2013 Garcia-Oliveira et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
BiolVlGCl Central Creative Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly cited. 



Garcia-Oliveira et al. BMC Plant Biology 2013, 13:134 
http://www.biomedcentral.eom/1 471-2229/1 3/1 34 



Page 2 of 13 



Background 

Aluminium (Al) toxicity is one of the major concerns for 
crop productivity in acidic soil, accounting for more 
than 50% of the global arable land [1]. Al is one of the 
highly abundant elements in the earth crust and under 
low pH conditions (acidic soils), it is solubilised in the 
soil in a toxic ionic form that inhibits root elongation. 
The primary response of plants to Al toxicity is the rapid 
inhibition of root growth, particularly in the root apex 
by blocking the process of cell division along with cell 
elongation, and subsequently inefficient absorption of 
nutrients and water from soil, resulting in the reduction 
of plant growth and overall productivity [2]. 

Plant species differ in the level of Al resistance because 
they evolved different mechanisms to overcome the se- 
lective pressure of Al toxicity imposed under acid soils 
[3]. These mechanisms can be broadly divided into two 
categories: Al resistance and Al tolerance. The Al resist- 
ance mechanism is often based on the Al exclusion from 
the root apex or the external detoxification of Al in the 
rhizosphere through the release of organic acid anions 
(oxalate, citrate and malate) from the roots limiting Al 
uptake [4]. On the other hand, the tolerance mechanism 
may involve the entrance of Al through the roots and its 
relocation (in the vacuole) or internal detoxification 
through Al chelating with organic acids (citrate and 
oxalate) [5]. These mechanisms have been validated at 
the molecular level, particularly with the functional 
characterization of the major genes such as ALMT 
(Aluminium-Activated Malate Transporter) and MATE 
(Multidrug and Toxic compound Exudation) in bread 
wheat (Triticum aestivum L.) [6] and sorghum (Sorghum 
bicolor L. Moench) [7], respectively. In addition, genes 
encoding transporters have been identified through mu- 
tational analysis, especially ABC transporters type (ALS1 
and ALS3) in Arabidopsis [8,9] and bacterial-type ABC 
transporters (STAR1 and STAR2) in rice (Oryza sativa 
L.) [10]. Recently, a plasma membrane localized trans- 
porter Nratl (Nramp aluminum transporter 1) was also 
found to be associated with Al tolerance particularly to 
trivalent form of Al in rice [11]. These genes are essen- 
tial for Al resistance, and seem to be a part of the path- 
way involved in the secondary level of protection 
through uptake and redistribution of Al to less sensitive 
tissues in plants, although their actual function in Al tol- 
erance mechanism in plant is still unclear. 

There is strong evidence supporting the important role 
of regulatory genes (transcriptional factors) in plant tol- 
erance to abiotic stresses [12-15]. Recently, the zinc fin- 
ger transcription factors, STOP1 (sensitive to proton 
rhizotoxicity) and ART1 (Al resistance transcription fac- 
tor) have been identified in Al sensitive mutants of 
Arabidopsis and rice, respectively [16,17]. Transcriptome 
analysis under Al stress revealed that 101 and 31 genes 



were down-regulated in the Arabidopsis stopl [18] and 
rice artl mutants [17], respectively. Interestingly, the 
major genes related to Al tolerance, particularly ALMT1 
and MATE1 were regulated by these transcription factors 
[18,19]. Recently, the contribution of ART1 locus to the 
variation for Al tolerance in rice has also been identified 
in QTL analysis [20]. The limited impact of single func- 
tional genes in plant stress tolerance has been associated 
with the polygenic nature of such traits. Thus, the identifi- 
cation and characterization of key regulatory genes that 
act as master regulators controlling entire response net- 
works would be the most promising and sustainable ap- 
proach to modify complex traits in plants as they 
coordinate the expression of many target genes [21], 

Wheat is one of the most important natural allopoly- 
ploid species, as it is not only directly or indirectly con- 
tributing in the food supply for nearly half of the worlds' 
population but also can serve as a model plant for other 
economically important polyploid crop species. It is con- 
sidered as one of the sensitive crop to Al stress among 
cereals. Bread wheat, with hexaploid nature comprised 
from three genomes (AABBDD) are organized in seven 
homoeologous groups, each homoeologous group has 
individual gene in triplicate form (from each of A, B and 
D genomes). Consequently, it is of great interest to re- 
veal how the expressions of homoeologues genes are 
regulated in hexaploid wheat because theoretically, all 
the three homoeologues of a gene are assumed to be 
uniformly expressed. In the previous decade, several 
studies have been conducted in order to identify the mo- 
lecular markers (random DNA markers) linked to Al 
tolerance through QTL mapping and genome-wide asso- 
ciation analyses [22-26]. So far only two candidate genes 
ALMT1 and MATE1 for Al tolerance in wheat have 
been identified and also mapped on chromosome 4DL 
using Chinese Spring deletion lines and 4BL through 
QTL mapping, respectively [22,27,28]. Recently, the 
three homoeologues of TaMATEl have been cloned 
[29], although, no information is available on the expres- 
sion of respective homoeologues of these candidate genes 
for Al tolerance in hexaploid wheat. Therefore, in order to 
improve Al tolerance in bread wheat, identification of 
three homoeologues of the candidate gene is a promising 
strategy that could be utilized to develop functional 
markers for genomic assisted breeding programme in 
wheat. Herein, we report on the identification, physical 
localization, and molecular characterization of a novel 
transcription factor TaSTOPl homoeologues genes in 
bread wheat for Al and proton (H + ) tolerance. 

Results 

Cloning and structure of TaSTOPl 

In order to clone the STOP1 in wheat, primers from 
wheat EST showing highest similarity with Arabidopsis 



Garcia-Oliveira et al. BMC Plant Biology 2013, 13:134 
http://www.biomedcentral.eom/1 471-2229/1 3/1 34 



Page 3 of 13 



thaliana STOP1 were used to amplify the TaSTOPl in 
bread wheat genotype Barbela 7/72/92. Further, 5' and 3' 
UTR ends of TaSTOPl were amplified using RACE as 
described in material and methods. TaSTOPl was ampli- 
fied in six different bread wheat genotypes and its multiple 
alignments suggested the mixed amplification from distinct 
wheat genomes (A, B, and D). Moreover, the comparison 
of the TaSTOPl cDNA sequence with genomic sequence 
revealed that TaSTOPl gene does not contain introns. 

The coding region of TaSTOPl (TaSTOPl-A) from 
genome A has 1533 bp length and is differentiated from 
genomes B (TaSTOPl-B) and D (TaSTOPl-D) by having 



a 6-bp deletion (1276-1281 positions). TaSTOPl-B and 
TaSTOPl -D are differentiated by several SNPs in the re- 
spective open reading frame (ORF) [Additional file 1]. 
TaSTOPl-D in bread wheat genotype Barbela 7/72/92 
has a full-length cDNA of 1970 bp, containing a coding 
region of 1539 bp that encodes a polypeptide of 512 
amino acids (Figure 1). TaSTOPl-D has a molecular 
weight of 55.8 kD, with a 5.64 isoelectric point (pi). 
InterProScan function domain analysis suggests that 
TaSTOPl belongs to the Cys2His2 zinc finger family 
protein. Subcellular prediction analysis indicated that 
TaSTOPl protein is localized in the nucleus. 
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Figure 1 Nucleotide sequence of TaSTOPl-D from Barbela 7/72/92 and its deduced protein (amino acid sequences). The nucleotide 
sequence was presented over the deduced amino acid sequence. Note: translation stop codon was noted by-. 
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Phylogenetic studies for ST0P1 gene among 
plant species 

Deduced amino acid sequence of TaSTOPl consists of 
510 (TaSTOPl-A) and 512 (TaSTOPl-B and TaSTOPl- 
D) amino acids in bread wheat and the phylogenetic 
relationship with other STOP-like proteins from 32 dif- 
ferent plant species was studied based on their full 
length sequences. Phylogenetic analysis of STOP like 
proteins clearly formed two distinct clusters which im- 
plied that STOP2 (AT5G22890) was completely distinct 
from STOP1 (AT1G34370) in Arabidopsis. Further- 
more, STOP1 and STOP2 protein from monocotyledons 
can be clearly distinguished from eudicots (Figure 2A). 
TaSTOPl homoeologues from bread wheat showed high 
similarity with STOP1 from Barley (Hordeum vulgare) 
and Brachypodium distachium with identities of 92 and 
87%, respectively. Noticeably, the phylogenetic relation- 
ship of TaSTOPl homoeologoues particularly TaSTOPl - 
A and TaSTOPl-D from bread wheat (AABBDD) 
displayed the maximum similarity with STOP1 like pro- 
teins from Triticum urartu (AA) and Aegilops tauschii 
(DD), respectively (Figure 2A). 

Sequence alignment of STOP1 protein from bread 
wheat with its homologous from cereals [Brachipodium, 
barley, rice, sorghum, maize (Zea mays) and foxmillet 
(Setaria italic)] and Arabidopsis, shows that TaSTOPl en- 
codes a putative Cys2His2 zinc finger protein containing 
four potential zinc finger domains. Three zinc finger do- 
mains (ZF1, ZF2, and ZF4) are predicted as the Cys2His2 
type, whereas ZF3 is predicted as the Cys2His-Cys or the 
Cys2His2 type (Figure 2B). Furthermore, sequence align- 
ment revealed that STOP-like proteins share highly con- 
served regions in the ZF domains. In addition, the STOP1 
protein also showed greater sequence conservation in the 
C-terminus than in the N-terminus among family mem- 
bers. The 6-bp deletion in TaSTOPl- A resulted in a loss 
of two amino acids [Pro (P) and Gin (Q)] at position 426 
without changing the reading frame [Additional file 2]. 
We also observed in Barbela 7/72/92, two SNPs at pos- 
ition 221 (A replaced by G) and 1186 (G replaced by A) 
which substituted an Asn (N) by a Ser (S) and Asp (D) by 
Asn (N) in TaSTOPl-B than in both TaSTOPl-A and 
TaSTOPl-D homoeologues, respectively [Additional file 
1]. Multiple alignment of STOP1 protein revealed that the 
SNP at position 1186 is located in the zinc finger domain 
(ZF4) of TaSTOPl homoeologues (Figure 2B). 

TaSTOPl localization on distinct wheat genomes 

For localization of TaSTOPl on homoeologous chro- 
mosomes, locus specific primer pairs were designed 
[Additional file 3] and TaSTOPl homoeologues were 
amplified in a series of nullitetrasomic lines of Chinese 
Spring wheat along with Chinese Spring as a positive 
control. On the basis of the presence or absence of 



PCR products visualized in agarose gels, we observed 
that TaSTOPl genes are located on homoeologous chro- 
mosomes 3 A, 3B and 3D, and named TaSTOPl- A, 
TaSTOPl-B and TaSTOPl-D, respectively (Figure 3A-C). 
Furthermore, ditelosomic lines of Chinese Spring for 
homoeologous group 3 chromosomes were also used to 
assign the TaSTOPl on chromosomal arms and confirmed 
that TaSTOPl-A, TaSTOPl-B and TaSTOPl-D genes are 
located on the long arm of homoeologous chromosomes 
3 A, 3B and 3D, respectively (Figure 3D). 

In order to further validate these results, physical 
chromosomal localization of TaSTOPl was performed 
using Tyramide Signal Amplification FISH (Tyr-FISH) 
technique in the root-tip metaphase chromosome spreads 
of bread wheat genotype Barbela 7/72/92. For this pur- 
pose, a probe was developed by PCR amplification of 
TaSTOPl from genomic DNA of Barbela 7/72/92. 
FISH with the TaSTOPl specific probe on metaphase 
chromosomes depicted hybridization signals of TaSTOPl 
on both chromatids of the long arms of three chromo- 
somes (Figure 4A) and simultaneous re-probing of 
chromosome preparations with GAA- (red) and pAsl 
(green) identified homoeologous group 3 with positive sig- 
nals on chromosomes 3AL, 3BL and 3DL (Figure 4B). 
Some background was also observed but it was distin- 
guished from true hybridization signals because they usu- 
ally appear as dots without any pattern and not in both 
chromatids of the same chromosome. Only those signals 
having the same size and appearance on the same position 
of both chromatids of one specific chromosome were ana- 
lyzed as true-positive hybridization signals. 

TaSTOPl transactivation activity 

The nuclear localization of TaSTOPl protein was pre- 
dicted by WOLF-PSORT programme [30]. TaSTOPl is 
a putative transcription factor of the Cys 2 His 2 -type Zinc 
fingers family. The three homoeologues of TaSTOPl 
showed similar genomic structure in bread wheat 
[Additional file 2]. Therefore, we only proceed to 
evaluate the transactivation potential of TaSTOPl lo- 
cated on genome A (TaSTOPl -A), To this aim, we 
performed a modified yeast one-hybrid assay [31] to 
evaluate the ability of the lexA-TaSTOPl-A fusion pro- 
tein to activate the lexA-driven expression of the lacZ 
gene in a heterologous yeast system, the results of which 
are illustrated in Figure 5. It is shown that lexA-TaSTOPl 
-A has the potential to transactivate lacZ expression either 
in the presence or absence of Al, further corroborating its 
potential role as a transcription factor. 

Expression profile of three TaSTOPl homoeologues in 
bread wheat 

To understand the homoeologue specific expression of 
TaSTOPl in bread wheat, we measured the temporal 
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Figure 2 Phylogenetic tree and multiple alignments of potential Zinc Fingers domain of STOP-like proteins. Phylogenetic tree based on 
amino acid sequences showing the relationship of TaSTOPI with other plant ST0P1 type proteins (A). The deduced amino acid sequences were 
aligned with CLUSTALW. Comparison of ZFs domain of TaSTOPI (TaSTOP 1 -A, TaSTOPI -B, TaSTOPI -D) with its homologous from Brachypodium 
(B.distachyon), Barley (H.vulgore), rice {O.sotivo), sorghum {S.bicolor), maize {Z.moys) foxmillet (S.itolico), and Arabidopsis {A.thalianal & A.tholiono2) 
(B). Horizontal bars indicate ZFs domain and asterisks indicate conserved motif of Cys2His2 or Cys2His2-Cys. 
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Figure 3 TaSTOPl mapping on homoeologous chromosomes and their arms in bread wheat. TaSTOPl mapping on wheat chromosomes 
from genome A (A), B (B) and D (C) using nullitetrasomic lines. M: Molecular-weight marker (100bp ladder), 1 to 7: Chinese Spring as control, 
N3AT3B, N3AT3D, N3BT3A, N3BT3D, N3DT3A, N3DT3B. Arms mapping of ToSTOPl-A, ToSTOPl-B and ToSTOPl-D with ditelosomic lines of 
homoeologous group 3 (D). 1, 4 and 7: Chinese Spring as a control; 2: Dt3AS; 3: Dt3AL; 5: Dt3BS; 6: Dt3BL; 8: Dt3DS and 9: Dt3DL 



expression of TaSTOPl-A, TaSTOPl-B and TaSTOPl-D 
by Real-time qPCR in the root and shoot tissues of two 
bread wheat genotypes [Barbela 7/72/92 (Al resistant) and 
Anahuac (Al sensitive)] grown under Al stress (74 uM Al). 
The 18S RNA expression was used as an internal control. 
We noticed a biased transcription of homoeologues of 
TaSTOPl gene in the root and shoot tissue of these di- 
verse bread wheat genotypes (Figure 6). The transcript 
levels of homoeologue TaSTOPl -A under control as well 
as Al treatment were significantly higher than those of 



TaSTOPl -B and TaSTOPl ~D in root and shoot tissues of 
both genotypes (depending upon genotypes, fold differ- 
ences in the transcript abundance were about 5-6 and 6- 
11 for homoeologue TaSTOPl-A versus TaSTOPl-B in 
root and shoot, respectively) whereas expression of 
homoeologue TaSTOPl -D was only two-fold higher 
than TaSTOPl-B in the root tissues of both genotypes 
(Figure 6 A and B). A slight induction (within the two 
hours of Al exposure) followed by return to basal levels 
of the homoeologue TaSTOPl -A transcripts was 




Figure 4 Physical localization of TaSTOPl on homoeologous chromosome in the root-tip spreads of bread wheat genotype Barbela 

7/72/92 using TYR-FISH technique. (A) Detection of positive hybridization signals of TaSTOPl on both chromatids of the long arms of three 

chromosomes (B) and identification of chromosomes after re-probing with GAA- (red) and pAs1 (green). Note: The GAA-satellite sequence 

identifies A and B genome chromosomes whereas the pAs1 identifies chromosomes from the D genome. Arrow indicates the position of 

TaSTOPl on respective chromosome (Scale bar = 10 microm). 
\ J 
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observed in the root tissues of Barbela 7/72/92 (Al re- 
sistant) whereas in the Anahuac (Al sensitive) under Al 
stress it was observed a quite stable transcript expres- 
sion (Figure 6A). The expression of TaSTOPl-B was 
also almost unaltered in the roots of both genotypes 
under Al stress. Similarly, TaSTOPl -D was also consti- 
tutively expressed under Al stress in the roots of both 
genotypes except a rapid (2 h) and significant repres- 
sion in the Al sensitive genotype Anahuac (Figure 6A). 
Nevertheless, a quite stable transcript expression of 
TaSTOPl homoeologues was observed in the shoot of both 
genotypes except homoeologue TaSTOPl -A (Figure 6B). 
Interestingly, and similarly to roots, the expression of 
TaSTOPl -A in the shoots of both genotypes was also 
slightly induced under Al stress, however, its up- 
regulation was noticed only after 6 h and return to its 
basal level (Figure 6B). It is worth noting that the tran- 
script levels of TaSTOPl -B of Al tolerant genotype 
Barbela 7/72/92 are lower in the shoot tissues, com- 
pared to the sensitive genotype Anahuac (Figure 6B). 

Moreover, to evaluate the transcript expression of 
TaSTOPl homoeologues under protons (H + ) stress, we 
also performed the relative quantification of homoeologue 
specific expression of TaSTOPl transcript in root tissues 
of both genotypes under different levels of pH (6.0, 5.0 
and 4.0). A slight gradual decrease in the transcript levels 
of TaSTOPl- A homoeologue was observed in the roots of 
Barbela 7/72/92 under protons (H + ) stress, but a quite 
stable expression was observed in genotype Anahuac 
(Figure 6C). Furthermore, considerable up-regulation of 
the homoeologue TaSTOPl ~B and slight repression of 
the homoeologue TaSTOPl -D transcript level was also 
observed under low pH (pH 4.0) in the roots of geno- 
type Anahuac and Barbela 7/72/92, respectively 
(Figure 6C). 

Discussion 

So far only ALMT1 and MATE genes have been de- 
scribed as responsible for most of the genotypic vari- 
ation for Al tolerance and are considered as major genes 
for Al resistance in bread wheat [32]. These genes have 
been identified in wheat [6,27] and also in the model 
plant Arabidopsis [19,33] and major cereals such as bar- 
ley [34], maize [35], sorghum [7] and rye [36,37]. The 
mechanisms underlying Al-tolerance in plant species 
have yet to be fully elucidated, as it seems that Al targets 
multiple cellular sites such as cell walls, plasma membranes, 
and cellular processes, like signal transduction path- 
ways and homeostasis mechanisms [2,5,38,39]. Recent 
evidence suggests that regulatory genes (transcription 
factors) play a key role in Al detoxification at different 
cellular levels [16,17]. 

In the present study, a novel transcription factor gene 
named TaSTOPl was cloned from an Al tolerant 



genotype Barbela 7/72/92 (Figure 1) which was derived 
from Portuguese bread wheat landrace Barbela through 
single seed descent method [40]. TaSTOPl belongs to a 
member of Cys2His2 zinc finger family proteins that con- 
tains four potential zinc finger domains and has highly 
conserved regions in the zinc finger (ZF) domains (Fig- 
ure 2B). Like other C2H2-type zinc finger proteins it con- 
tains more than one zinc finger motif with highly 
conserved amino acid sequence for DNA binding 
[16,41]. Phylogenetic analysis clearly differentiated the 
STOP1 like proteins from monocots in a group (Fig- 
ure 2A) and also suggested that genome A and B are 
more distant from genome D [Additional file 2], as gen- 
ome D was incorporated recently in comparison to the A 
and B genomes in bread wheat [42]. In addition, our 
phylogenetic analysis results also support previous find- 
ings at molecular level that the A and D genomes in bread 
wheat were derived from T urartu and A. tauschii, re- 
spectively [Additional file 2] [42]. 

In past, the nullitetrasomic and ditelosomic lines of 
Chinese Spring wheat have been extensively used in clas- 
sical and molecular genetic studies for the identification 
of loci associated with numerous traits including Al 
tolerance [25,43]. Thus, gene mapping has practical im- 
plications in plant genetics and breeding, and the deter- 
mination of its physical localization is even more precise 
as it enables us to confirm the exact position of a gene 
on a chromosome. PCR based mapping of TaSTOPl 
using nullitetrasomic and ditellosomic lines of Chinese 
Spring wheat revealed that TaSTOPl is located on the 
long arm of wheat homoeologous group 3 chromosomes 
(Figure 3). 

In situ hybridization is not only the most direct 
method for physical localization of genes in chromo- 
somes, but also it is a valuable tool for the identification 
of copy numbers of a gene in species having a highly 
complex genome such as the bread wheat. Due to its 
hexaploid nature, most genes can be found in triplicate 
with one copy on each genome [44-46]. Tyramide signal 
amplification FISH (Tyr-FISH) technique has been suc- 
cessfully used for the identification of low copy number 
DNA sequences in wheat [47]. However, for the success- 
ful localization of a gene using the FISH technique, the 
minimum length of the probe size seems to be crucial. 
Recently, RD50 gene has been localized in bread wheat 
with a probe size of 2 kb [47]. Physical localization of 
TaSTOPl by Tyr-FISH technique not only confirmed 
the results obtained from PCR based mapping, but also 
indicated that TaSTOPl could be a single copy gene 
localized on the long arm of homoeologous group 3 
chromosomes (Figure 4). In bread wheat, physical 
localization of single copy gene Glu-1 and RDS0 with 
FISH technique has been demonstrated on long and 
short arms of the homoeologous group 1 chromosomes, 
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respectively [47,48]. In the present work, we were able 
for the first time to localize a gene with a probe size 
smaller than 2.0 kb demonstrating that the FISH technique 
can be used to simultaneously anchor homoeologous chro- 
mosomes with 1.5 kb probes even in bread wheat. 

Subcellular prediction of TaSTOPl protein in the nu- 
cleus is in agreement with data previously reported [18]. 
Zinc finger motifs are thought to recognize and bind to 
target DNA sequences, but they are not required for tran- 
scriptional activity [49,50]. Our results clearly exhibited 
the transactivation potential of TaSTOPl -A at least in 
yeast (Figure 5). In present investigation, we observed 
highly similar genomic structure of TaSTOPl genes 
[Additional files 1 and 2], therefore, it seems that 
TaSTOPl transactivation function is constitutive and 
may not depend on the presence/absence of Al. The 
SNPs observed among TaSTOPl homoeologues showed 
minor changes in respective amino acids of putative 
proteins which could alter the secondary structures by 
influencing the folding of these proteins (Additional file 1). 
Therefore, the three homoeologues of the same locus in 
bread wheat share high sequence similarity that could illus- 
trate the flexibility of a polyploidy species in which due to 
the mutation in one homoeologue may be compensated 
for by the homoeologue [29]. 

Allopolyploidy plays an important role in plant evolu- 
tion that arises with the merging of two or more ge- 
nomes into single nucleus which may contribute either 
equally or disproportionately. However, recent molecular 
findings confirmed the asymmetric genomic expression 
pattern in natural and synthetic allopolyploid plant spe- 
cies, but the predominant transcript expression of one 
genome over the other genome(s) vary from gene to 
gene [51]. In hexaploid wheat, a uniform level of expres- 
sion for all the three homoeologues has been reported 
for approximately 20% of the unigene loci [52], whereas 
20-29% genes did not express at least one homoeoallele 
[52,53]. Thus, the relative contribution of the three 
homoeologues of TaSTOPl gene at transcript expression 
level was determined under Al and proton (H + ) stresses 
in two bread wheat genotypes showing contrasting be- 
haviour for Al toxicity. Transcript expression profiling of 
TaSTOPl homoeologues in root and shoot tissues iden- 
tified the predominance of TaSTOPl -A homoeologue 
followed by TaSTOPl -D over TaSTOPl -B in root and 
only predominance of TaSTOPl- A in shoot tissues of both 
genotypes under control and stress (Al and pH) condi- 
tions (Figure 6). Similarly, among the three homoeologues, 
the higher expression of one homoeologue of a MAD box 
transcription factor and Spa gene have also been observed 
in bread wheat [54,55]. Although, the presence of cis ele- 
ments within the 5' untranslated region of a gene is un- 
usual, it is not aberrant or abnormal [56]. Surprisingly, 
sequence analysis of 5' UTR of TaSTOPl genes 



differentiated the TaSTOPl -A homoeologue from 
TaSTOPl ~B and TaSTOPl -D due to the presence of a 
pyrimidine-rich stretch and absence of light responsive 
element [Additional file 4]. In tomato, the deletion of 
the 5' UTR containing pyrimidine-rich stretch from the 
HMG2 promoter reduced the level of HMG2 gene ex- 
pression by a factor of 10 [56]. Therefore, among the 
three homogeologues in bread wheat, the up- or down- 
regulation of the expression of specific homoeologue 
could be the result from either the dominancy of ances- 
tral diploid donor parent or early polyploidization-cis- 
regulatory variation. 

Interestingly, the time-dependent Al-responsive ex- 
pression of TaSTOPl homoeologues observed in the 
root tissues of two bread wheat showing contrasting 
phenotypes for Al toxicity suggests a putative role for 
TaSTOPl in Al resistance (Figure 6 A and B). Al respon- 
sive STOP1 expression has also been reported in 
Arabidopsis [16], alfalfa [57] and common bean [58]. 
Similarly to Arabidopsis, the transcript expression of 
homoeologues of TaSTOPl in the roots of diverse bread 
wheat genotypes was also modulated in response to 
proton stress [16]. It is noticeable that genotype Barbela 
7/72/92 is highly resistant to Al toxicity compared with 
Anahuac, but in the absence of Al under low pH we 
did not observe significant differences between these 
genotypes for root growth [Additional file 5]. It is worth 
noticing that the three homoeologues of TaSTOPl have 
similar genomic structures, but showed a different tran- 
script expression in response to Al and proton (H + ) stress. 
These results may reveal that the homoeologues of 
TaSTOPl are differentially contributing to Al or proton 
(H + ) tolerance in bread wheat. Further work is in progress 
in order to definitely establish if these genes play a signifi- 
cant role in aluminium tolerance. 

Conclusions 

Bread wheat has wide genotypic variation for Al resist- 
ance [23,40] and the role of several chromosomes such 
as chromosome arms 2DL, 3DL, 4BL, 4DL, 6AL, 7AS 
and chromosome 7D, in Al tolerance has been revealed 
in classical genetic studies through chromosome ma- 
nipulation in wheat [43]. Furthermore, several QTL as- 
sociated with Al tolerance in bread wheat have also been 
reported by many researchers [22,24-26]. Contrarily to 
Arabidopsis and rice, due to the paucity of a priori can- 
didate genes in wheat only two major QTL located on 
chromosome 4DL and 4BL have yet been elucidated 
at molecular level, showing the co-segregation with 
candidate genes TaALMTl homoeologue (4DL) and 
TaMATEl homoeologue (4BL), respectively [22,27]. 
The classical studies using chromosomal manipula- 
tion as well as recent QTL mapping and genome-wide 
association analysis have also detected loci associated 
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with Al resistance on homoeologous group 3 chromo- 
somes (3A, 3B and 3D) in bread wheat [23-25,43]. Re- 
cently, the role of a zinc finger transcription factor 
ART1 identified through mutational analysis in rice 
has also been shown in natural variation of Al toler- 
ance in rice, earlier which was suggested that it was 
not involved in Al tolerance [20,32]. 

In the present investigation, we cloned and character- 
ized the novel candidate genes TaSTOPl in bread wheat. 
The homoeologues of TaSTOPl exhibited similar gen- 
omic structures, but showed biased transcript expression 
and different response to Al and proton (H + ) toxicity. 
Furthermore, TaSTOPl -B homoeologue from Al tolerant 
genotypes Barbela 7/72/92 and Viloso Mole not only 
showed highest similarity but also contain the same SNP 
located in ZF4 domain than Al sensitive genotypes 
Anahuac, Chinese Spring and Saloio [Additional file 2] 
[40]. Finally, in order to correlate with Al tolerance, it 
would be very interesting either to functionally characterize 
or further verify the role of TaSTOPl because gene(s) 
underlying the QTL on homoeologous group 3 chromo- 
somes particularly 3BL has not been so far identified in this 
important cereal. 

Methods 

Plant material and growth condition 

The seedlings of selected bread wheat genotypes Barbela 
7/72/92, Anahuac, Chinese Spring, Ruivo, Viloso Mole 
and Saloio classified in relation to Al tolerance [40] were 
grown in hydroponic solution (0.4 mM CaCl 2 , 0.65 mM 
KN0 3 , 0.25 mM MgCl 2 .6H 2 0, 0.1 mM (NH 4 ) 2 S0 4 and 
0.04 mM NH4NO3) and kept in controlled growth 
chamber under 14 h/26°C day and a 10 h/22°C night re- 
gime, with a light intensity of 150 \imo\ photons m "V 1 
and a relative humidity of 65%. 

Al stress assay 

Four days old seedlings of Al resistant bread wheat 
genotype Barbela 7/72/92 and Al susceptible genotype 
Anahuac growing in hydroponic solution with pH 4.0 
were shifted in fresh nutritive solution having 74 \iM Al 
whereas control seedling were raised in only fresh 
hydroponic solution. For root re-growth measurement, 
after 24 h, roots were immersed in 0.1% eriochrome cya- 
nine R (Sigma, Germany) dye solution for 10 minutes and 
prior to measurement; seedlings were allowed to grow in 
fresh nutritive solution for 48 h. For visual detection of Al 
accumulation in roots, hematoxylin assay was performed 
[59] and root samples were observed under microscope. 

DNA and RNA extraction 

Total genomic DNA was isolated from young leaf sam- 
ples using the DNeasy Plant Mini Kit (Qiagen, Germany) 
and total RNA was extracted using Trizol method 



followed by purification using PureLink™ RNA Mini Kit 
(Ambion, Invitrogen, USA). The first-strand cDNA was 
synthesized in a final volume of 20 \A reaction containing: 
1 ug RNA, 2 ul lOx RT buffer, 0.8 [A of 25x dNTP mix 
(100 mM), 2 (il of lOx RT random primers and 1 \A of 
Multiscribe reverse transcriptase (Applied Biosystems, 
USA). 

Gene structure and cloning of full-length cDNA 

The Arabidopsis thaliana AtSTOPl (AT1G34370) pro- 
tein sequence was retrieved from TAIR database (http:// 
www.arabidopsis.org) and was used as the query se- 
quence to search the GenBank wheat expressed se- 
quence tag (EST) database. The homologous sequence 
of STOP1 in wheat was obtained via tblastn in NCBI 
(http://www.ncbi.nlm.nih.gov). Homologous wheat EST 
clones showing highest similarity with Arabidopsis STOP1 
were retrieved. TaSTOPl (wheat EST) specific primers 
were designed to perform the sequential analysis in bread 
wheat genotype Barbela 7/72/92 [Additional file 3]. Total 
RNA from the seedlings of Barbela 7/72/92 and Anahuac 
was employed to obtain the full-length cDNA of TaSTOPl 
including 5' UTR and 3' UTR using Rapid Amplification 
of cDNA Ends (RACE) (SMARTer™ RACE cDNA Ampli- 
fication Kit, clontech, USA). The genomic sequence of 
TaSTOPl was obtained in six different bread wheat geno- 
types, namely Barbela 7/72/92, Anahuac, Ruivo, Viloso 
Mole, Saloio and Chinese Spring [Additional file 6]. 

Molecular analysis and construction of phylogenetic tree 

The Pfam (http://www.sanger.ac.uk/Software/Pfam/search, 
shtml) software was used to identify potential domains and 
WOLF-PSORT (http://www.psort.org) was used to predict 
the intracellular localization of the TaSTOPl protein. Fifty 
two STOP1 like protein sequences were collected from 32 
different plant species including TaSTOPl homoeologues 
from bread wheat genotype Barbela 7/72/92 using the Basic 
Local Alignment Search Tool (BLAST) programme from 
NCBI (http://blast.ncbi.nlm.nih.gov/Blast.cgi) [Additional 
file 7]. Phylogenetic tree was constructed by MEGA4 soft- 
ware using neighbor- joining method [60]. 

Mapping of TaSTOPl gene in wheat genome 

All the nullitetrasomic and ditelosomic lines of Chinese 
Spring wheat except for the homoeologous group 2, 4 and 
6 chromosomes used in present investigation were avail- 
able in our lab, whereas remaining stocks were kindly pro- 
vide by Prof. Adam J. Lukaszewski at Department of 
Botany and Plant Sciences, University of California, USA. 
On the basis of TaSTOPl alignment, different pairs of spe- 
cific primers for each TaSTOPl homoeologue were 
designed [Additional file 3]. TaSTOPl homoeologues 
(TaSTOPl-A, TaSTOPl-B and TaSTOPl-D) were ampli- 
fied in the nullitetrasomic (N3AT3B, N3AT3D, N3BT3A, 
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N3BT3D, N3DT3A and N3DT3B) and ditelosomic [Dt 
3AL and 3AS (long and short arm), Dt 3BL and 3BS, Dt 
3DL and 3DS] lines of Chinese Spring wheat. TaSTOPl 
homoeologues were assigned to the chromosomal arms 
based on the presence/absence of the PCR amplification 
products. 

Physical localization of TaSTOPl was performed by in 
situ hybridization. Probe was prepared by PCR amplifi- 
cation of the 1.5 kb genomic region of TaSTOPl from 
Barbela 7/72/92. Chromosome spreads from the root 
tips of Barbela 7/72/92 seedlings at mitotic metaphase, 
probe labelling and hybridization mixture were carried 
as described by Prieto et al [61]. Probe labelling was 
confirmed by dot-blot and detection of hybridization 
signals was carried out using the Tyramide Signal Amp- 
lification Kit (TSA™, PerkinElmer Life and Analytical 
Sciences, Inc., Waltham, MA, USA). To identify wheat 
chromosomes with positive signals, samples were re- 
hybridized using simultaneously the pASl repetitive se- 
quence and GAA-satellite sequence as probes [62,63]. 
The GAA-satellite sequence identifies all the A and B 
wheat chromosomes [62] whereas the pAsl identifies 
chromosomes from the D genome [63]. Individual 
slides were observed under a Nikon Eclipse 80i micro- 
scope (Nikon Instruments Europe BV, UK). Images 
were captured with a Nikon CCD camera using the ap- 
propriate Nikon 3.0 software and processed with 
Photoshop 4.0 software (Adobe Systems Inc., San Jose, 
California, USA). 

Transactivation assay 

TaSTOPl was amplified by PCR from Barbela 7/72/92 
using oligonucleotide primers containing restriction site 
for Kpn I [Additional file 3]. The Kpn I digested 
TaSTOPl fragment was cloned in frame with the lex A 
gene in the plasmid YCp91 [31]. Transformation of the 
ligated products was performed in the Escherichia coli 
strain XLl-Blue recAl endAl gyrA96 thi-hsdR17 supE44 
relAl lac [F'proAB lacPZDMIS TnlO (Tef)] (Stratagene, 
USA) and the positive clones were sequenced. The re- 
combinant plasmid YCp91-TaSTOPl-A was transformed 
in the yeast strain EGY48, which harbours the plasmid 
pSH18.84 encoding the lacZ gene under the control of a 
promoter containing LexA - responsive cis elements. 
Transformants were grown in selective media and |3- 
galactosidase measurements were performed as previ- 
ously described [64,65]. 

Analysis of TaSTOPl expression level 

For Al assay, seedlings of two bread wheat genotypes 
Barbela 7/72/92 and Anahuac showing contrasting re- 
sponse to Al toxicity, were raised in hydroponic nutrient 
solution with pH 4.0 during four days and further, were 
shifted to fresh nutrient solution with (74 \iM Al in the 



form of A1C1 3 -6H 2 0; stress treatment) or without Al 
(control treatment). Both root and shoot tissues were 
collected separately after treatment at specific time 
points (0, 2, 6, 12 and 24 h) from control and Al stress 
imposed seedlings. For pH assay, four days old seedlings 
of both genotypes grown in hydroponic nutritive solu- 
tion with pH 6.0 were transferred in fresh nutritive solu- 
tions having different levels of pH (pH 6.0, pH 5.0 and 
pH 4.0). Root samples of both genotypes from each 
treatment were collected after 0 h and 6 h exposure to 
different levels of pH. 

Three biological replicates of each sample were pre- 
pared and duplicate quantitative assays were performed 
for each cDNA sample. TaSTOPl homoeologues gene 
expression pattern was determined using the SYBR 
Premix Ex Taq (Takara, Japan) and the ABI 7500 Real- 
Time FAST PCR System (Applied Biosystems, USA). 
The 2 AAC T method [66] was used to quantify the relative 
expression levels of TaSTOPl homoeologues in compari- 
son to 18SRNA endogenous control. 

Accession numbers 

Sequence data from this article can be found in the 
GenBank database under the following accession num- 
bers: Barbela 7/72/92 (TaSTOPl-A: GenBank number 
KF034793; TaSTOPl -B: GenBank number KF034794; 
TaSTOPl-D: GenBank number KF034795), Anahuac 
(TaSTOPl-A: GenBank number KF034796; TaSTOPl-B: 
GenBank number KF034797; TaSTOPl-D: GenBank num- 
ber KF034798), Viloso Mole (TaSTOPl-A: GenBank num- 
ber KF034801; TaSTOPl-B: GenBank number KF034802; 
TaSTOPl-D: GenBank number KF034803), Saloio 
(TaSTOPl-A: GenBank number KF034804; TaSTOPl-B: 
GenBank number KF034805; TaSTOPl-D: GenBank num- 
ber KF034806), Chinese Spring (TaSTOPl-A: GenBank 
number KF034799; TaSTOPl-B: GenBank number 
KF034800) and Ruivo (TaSTOPl-D: GenBank number 
KF034807). 
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Additional file 1: Multiple alignments of the homoeologues of 
TaSTOPl in bread wheat genotype Barbela 7/72/92. 

Additional file 2: Phylogenetic tree and multiple alignments of 
TaSTOPl homoeologues in different species of wheat including 
some bread wheat genotypes. 

Additional file 3: Detail of primers used in present investigation. 

Additional file 4: Alignment of homoeologues of TaSTOPl 
5UTR region. 

Additional file 5: Physiological characterization of Barbela 7/72/92 
and Anahuac roots under Al stress. 

Additional file 6: Nucleotide and protein sequence of TaSTOPl in 
different bread wheat genotypes. 

Additional file 7: Detail of STOP like proteins in different plant 
species used for phylogenetic analysis. 
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